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Abstract 

i-£h ■ The radiation reaction for a point-like charge coupled to a massive scalar field 

is considered. The retarded Green's function associated with the Klein-Gordon 
wave equation has support not only on the future light cone of the emission point 
(direct part), but extends inside the light cone as well (tail part). Dirac's scheme of 
decomposition of the retarded electromagnetic field into the "mean of the advanced 
and retarded field" and the "radiation" field is adapted to theories where Green's 
function consists of the direct and the tail parts. The Harish-Chandra equation of 
motion of radiating scalar charge under the influence of an external force is obtained. 
This equation includes effect of particle's own field. The self force produces a time- 
changing inertial mass. 

1 Introduction 

There is an extensive literature devoted to the regularization problem in curved spacetime. 
De Witt and Brehme [1] derived the self-force acting upon a point electric charge in fixed 
background gravitational field. (Their expression for radiation reaction was later corrected 
by Hobbs [2].) The principle of equivalence in general relativity implies that a particle of 
infinitesimal mass and size moves along a geodesic. The particle of a small but finite mass 
perturbs the spacetime geometry. Mino, Sasaki and Tanaka [3] included the interaction of 
the particle with this perturbation which changes the worldline. In [4] Quinn has obtained 
an expression for the self-force on a point-like particle coupled to a massless scalar field. 
There is a long series of papers devoted to the motion of a point electric charge, a point 
scalar charge, and a point mass in black hole spacetimes, where the effects of radiation 
reaction are taken into account (see references in review [5]). 

The computation of effect of particle's own field is not a trivial matter, since the 
Green's function associated with the wave operator has support within the light cone. 
This is because in curved spacetime massless waves propagate not just at a speed of light, 
but also at all speeds smaller than or equal to the speed of light. The particle may "fill" its 
own field, which acts on it just like an external one. The equation of motion require one 
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to identify that portion of the retarded field at each point of the world line which arises 
from source contributions interior to the light cone. This part of field is often called the 
"tail term". The self force on a particle then consists of two parts: this comes from the 
direct part of the Green's function and depends on the current state of particle's motion 
and that comes from the tail part and depends not only the current state of the particle, 
but also on its past history. It leads to the non-local (integro-differential) equations of 
motion. 

Detweiler and Whiting proposed [6] a consistent decomposition of the retarded Green's 
function into singular and radiative parts. It obeys the spirit of Dirac's scheme of split- 
ting of electromagnetic potential of a point-like charged particle arbitrarily moving in 
flat spacetime. Dirac [7] decomposed the retarded Lienard-Wiechert potential A vct into 
two parts: (i) one-half of the retarded plus one-half of the advanced potentials which is 
inhomogeneous solution of the wave equation DA a = —Anxj a whose source term is infinite 
on the world line. A s = l/2(A mt + A adv ) is just singular as the retarded potential in the 
immediate vicinity of the particle's world line. The superscript "S" stands for "singular" 
as well as "symmetric", (ii) combination A R = l/2(A ret — A adv ) of one-half of the retarded 
minus one-half of the advanced potentials which satisfies the homogeneous wave equation. 
This well behaved potential can be thought as a free radiation field. The superscript "R" 
stands for "radiative" as well as "regular" . 

The radiative Green's function implicitly used by Dirac in flat spacetime is 

G™ d (x,y) = G* ct (x,y)-G s ^(x,y) (1.1) 

where 

G s ^{x,y) = \ [G ret (x,y) + G adv (x,y)] . (1.2) 

The causal structure of the Green's function is richer in curved spacetime. Due to 
contributions of the interior of the light cones, the retarded potential depends on the 
particle's history prior to the retarded instant T rct (x) while the advanced one is generated 
by portion of particle's world line ( after the advanced instant (The retarded and 

the advanced moments label the points on ( related with arbitrary field point x G M 4 
by null rays.) The combination of half- retarded potential minus half-advanced one could 
satisfy the homogeneous wave equation. Moreover, it would be smooth on the world 
line. But a self-force constructed from this radiative potential will be highly non-causal. 
It will be depend on particle's entire history, both past (through the retarded Green's 
function) and future (through the advanced Green's function). The Dirac's scheme (1.1) 
for decomposition cannot be adopted without modification in curved spacetime. The 
modification is performed in Ref.[6]. 

Detweiler and Whiting start with a Hadamard construction of a symmetric scalar field 
Green's function 

G s ^(x, y) = i- [u(x, y)S(a) - v(x, y)9(-a)} , (1.3) 

where u(x,y) and v(x,y) are smooth functions of the base point y and the field point x 
and a is half of the square of the distance measured along the geodesic from x to y. The 
step function 9 (—a) means that this part of G sym (x,y) has support within both the past 
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null cone and the future null cone of x. Having coupled it with S— shaped distribution of 
scalar charge moving along a geodesic we obtain the tail part of the field ip as follows: 



C0*0 = - 



drv(x, z(t)) + 



'(*) 



drv(x, z{t)) 



;i.4) 



To remove the noncausality, authors add to singular Green's function (1.3) biscalar v(x, y) 
being solution of homogeneous wave equation. A new symmetric (singular) Green's func- 
tion 



G s (x, V) = ^ [u(x, y)S(cr) - v(x, y)9(<j)} , 



[1.5) 



has no support within the null cone. Corresponding tail part of the scalar field depends 
on particle's history during the interval [r rct (x), r adv (x)]: 



'(*) 



drv(x, z(t)) 



T ret (x) 



(1.6) 



It is inhomogeneous solution of wave equation. Since it is singular just as the retarded 
tail field ^[^(x), the authors define the radiative Green's function as follows 



G K (x,y) = G rct (x,y)-G s (x,y) 



(1.7) 



(cf. eq.(l.l)). Corresponding radiation field 



^ R (x) = 



qu{x,z{r)) 
2& 



T adv (l) 



T rct (x) 



/ drv (x, z{t)) + - / dTv(x,z(r)) 

J -oo A Jr Ict (x) 



T adv (» 
(X) 



;i.8) 



is well behaved solution of homogeneous wave equation. In the coincidence limit, where 
field point x approaches to the world line at point z(t), the interval [r rct (x) , r adv (x)} 
shrinks to zero. Due to taking of this limit the potential (1.8) is generated by the portion 
of the world line that corresponds to the causal interval ] — oo, r]. The last integral in eq. 
(1.8) gives no contribution to a self-force because it cancels the ill defined part of gradient 
A,,^] coming from implicit dependence of r rct upon x. 

Following this scheme, Detweiler and Whiting recovered the results [1, 3, 4] for elec- 
tromagnetic, scalar, and gravitational fields. 

It is obvious that the physically relevant solution of the wave equation is the retarded 
potential. Teitelboim [8] derived the electromagnetic self-force in flat spacetime within the 
framework of retarded causality. The author substituted the retarded Lienard-Wiechert 
field in the Maxwell energy-momentum tensor density and calculated the flow of energy- 
momentum which flows across a space-like surface. Minkowski space was parameterized 
by four curvilinear coordinates. The first, proper time, labels points of emission placed on 
(, the second one determines the surface (e.g., a tilted hyperplane which is orthogonal to 
particle's 4- velocity at fixed instant of observation). Having integrated the stress-energy 
tensor over two angular variables that distinguish points on the surface, Teitelboim found 
the flow of energy-momentum mentioned above. The resulted expression depends on the 



3 



particle's individual characteristics (on its mass, its charge, its velocity and acceleration). 
In fact, the surface integration is equivalent to taking of coincidence limit in Dirac's 
scheme. Abraham-Lorentz-Dirac expression for electromagnetic self-force is obtained in 
[8] via consideration of energy-momentum conservation. 

Teitelboim [8] demonstrates that each of two terms which constitute Abraham radiation- 
reaction four- vector originates from the specific part of Maxwell energy-momentum tensor 
density. The author decomposes the stress-energy tensor into two divergent-free compo- 
nents: T = T bn d + T rad . Surface integration of the bound part T hnd results the flow which 
never gets far from the point-like source. It is permanently "attached" to the charge and 
is carried along with it. A charged particle cannot be separated from its bound elec- 
tromagnetic "cloud" which has its own 4-momentum and angular momentum. "Bare" 
charge and electromagnetic "cloud" constitute new entity: dressed charged particle. T bnd 
contributes into particle's inertia: 4-momentum of dressed charge e contains, apart from 
usual velocity term, the extra term whose time derivative is exactly the negative of the 
Schott term: 

2e 2 

= m/ a". (1.9) 

3 

Calculation of flux of the radiative component T rad yields the Larmor relativistic rate 
of radiated energy-momentum. This part of energy-momentum detaches itself from the 
charge and leads an independent existence. 

In this paper we tend Teitelboim's approach to allow contributions from interior of 
the light cone. For the clearest demonstration of the impact of our analysis, we refer 
to a point-like particle of mass m and charge e coupled to electromagnetic field in flat 
spacetime of three dimensions [9, 10]. In 2 + 1 electrodynamics the retarded Green's 
function associated with D'Alembert operator □ = rf a ^d a d^ is supported within the light 
cone [11, 12]: 

G M^°7:- |x - yi> (i.io) 

6(x° — y° — |x — y|) is the light cone step function defined to be one if x° — y° > |x — y| 
and defined to be zero otherwise. Synge's world function in flat space-time is numerically 
equal to half the squared distance between x and y: [5] 

a(x,y) = ^(x a -y a )(x^-y^). (1.11) 

The analysis of the simplest model with tails will give more deep understanding of De- 
tweiler and Whiting scheme of decomposition. 

The paper is organized as follows. In Section 2 and 3, we define the scalar poten- 
tials and scalar field strengths. In Section 4 we split the energy-momentum and angular 
momentum carried by massive scalar field into bound and radiative parts. Extracting 
of radiated portions of Noether quantities is not a trivial matter, since the massive field 
holds energy and momentum near the source. In Section 5, we derive equation of motion 
of radiating scalar pole via analysis of energy-momentum and angular momentum bal- 
ance equations. It is of great importance that conservation laws yield the Harish- Chandra 
equation of motion [13]. In Section 6, we discuss the result and its implications. 
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2 Green's function and potentials 



The dynamics of a point-like charge coupled to massive scalar field is governed by the 
action [14, 5] 

total — /part + ^int + Afield- (2.1) 

Here 

/field = ~^j d ' x + A: V) (2.2) 

is an action functional for a massive scalar field ip in flat spacetime with metric tensor 
i] al3 = diag(— 1, 1, 1, 1). The mass parameter k is a constant with the dimension of 
reciprocal length. Its physical sense will be discussed below. The integration is performed 
over all the spacetime. The particle action is 



/part = -m J drV-i 2 (2.3) 

where mo is the bare mass of the particle which moves on a world line ( G M 4 described by 
relations z a {r) which give the particle's coordinates as functions of proper time; z a {r) = 
dz a (r)/dT. Finally, the interaction term is given by 

/mt = 9 J drv^y^) (2.4) 

where g is scalar charge carried by a point-like particle. 

Variation on field variable ip of action (2.1) yields the Klein-Gordon wave equation 

{n-k 2 ) V (x) = -4np(x), (2.5) 

where □ = rf^d a d$ is the D'Alembert operator. We consider a scalar field (p(x) with a 
point particle source 

/+oo 
dr5\x-z(r)), (2.6) 
-oo 

where g is coupling constant and S A (x — z(t)) is a four-dimensional Dirac's distribution 
localized on the world line: charge's density is zero everywhere, except at the particle's 
position where it is infinite. 

The relevant wave equation for the Green's function G(x, y) is 

{U-kl)G{x-y) = -5\x-y). (2.7) 
Its solution is the symmetric Green's function [5, 14] 

1 



G""(*-») lrr 



5(a) - -^=.h(k ^2^)e(-o) 



(2.8) 



where J\ is the first order Bessel's function whose argument contains Synge's world func- 
tion (1.11). 
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Convoluting the retarded Green's function 



9(x° - y°)G sym (x - y) 



(2.9) 




Mi(WT 2 - W) 



VT 2 - R 2 



where T = x° — y° and R — |x — y| with charge density (2.6), we construct the massive 
scalar field. The delta-function in eq.(2.9) results direct term which is generated by a 
single event in space-time: the intersection of the world line and x's past light cone. The 
Heaviside stepfunction extends the support to the segment of ( that corresponds to the 
interval r e] — oo, T TCt (x)]. Summing up the direct term and the tail term, we construct 
the retarded potential produced by an arbitrarily moving point-like source coupled to 
massive scalar field [14, 13, 15]: 



The first term is referred to the retarded point z rct associated with x; r is the retarded 
distance between z vct and x. K = x — z{r) is four-vector pointing from z{r) G ( to a 
point x where the potential ip is observed (see Figure 1). 



Figure 1: Direct part of the retarded potential at field point x is generated by the inter- 
section z(T rct ) of the world line ( and past light cone with vertex at x. Tail part of ip ret (x) 
depends on the segment of ( that corresponds to time interval ] — oo, r rct (x)]. The tail 
term defines contribution from cone's interior. The vector K is a vector pointing from 
the point of emission z(t) G ( to a field point x. 

Convolution of the advanced Green's function 




(2.10) 



— oo 




c 



'adv 



(x-y) 



6(-x° + y°)G sym (x-y) 



(2.11) 




k Ji(k VT 2 - R?) 



VT 2 - R 2 
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coupled with the 5- like density (2.6) is analogous to the manipulations with its retarded 
counterpart. The advanced potential 

r dT Mi^Ei ( , 12) 

^adv J yf-{K-K) 

T* dv (x) 

is generated by the point charge during its entire future history following the advanced 
time associated with x (see Figure 2). Particle's characteristics in the direct term are 
referred to the instant r adv (x). 




Figure 2: The advanced Green's function (2.11) possesses support on the surface of the 
future light cone of x (direct part) as well as within this cone (tail part). The portion 
of the world line after r adv (x) produces the advanced potential (2.12). Advanced instant 
labels point at which particle's world line punctures future light cone with vertex at x. 

The simplest scalar field is generated by an unmoved source placed at the coordinate 
origin [17, 18]: 

exp(-/c r) {n , q\ 

<Pyu{x) = g . (2.13) 

r 

Here r = a/ (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 is the distance to the charge. It is the well-known Yukawa 
potential which plays an important role in the theory of nucleon-nucleon interactions. 



3 Scalar field strengths 

Let two particles interact through a massive scalar field. Static charge g\ placed at the 
coordinate origin generates the Yukawa field with components 
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= -gi-^-(l + kor)ni, 

where rii = Xi/r is unit direction vector and r is the distance to the charge. This charge 
exerts another one, say #2, with the force 

g -fcori2 

F12 = -9i92 — — (1 + k r) n 12 , (3.2) 

where rii2 is unit vector codirectional with the radius vector ri2 drawing from particle 1 
to particle 2. The force is attractive if g\ and (?2 are of like sign. The charges can be 
interpreted as nucleons which attract each other and are bound in a nucleus. The force 
(3.2) elucidates also why a neutron repels an antineutron. 

In general, the action (2.1) can be considered as a classical version of the Yukawa 
model for strongly interacting nucleons. It is assumed that the neutrons and protons are 
joint together by a massive pseudoscalar pion field. Parameter ko associates with the rest 
mass of a massive scalar field particle mediating the interaction. The parameter acts as 
a cutoff: the heavier is the field particle, the shorter is range of Yukawa force. Yukawa 
force is essential at a distance about l/k . 

Let us derive the field produced by an arbitrarily moving scalar charge. Scalar field 
strengths are given by the gradient of the potential (2.10). Differentiation of the direct 
term yields 



[1 + (K ■ a)]K 



ii 



^"3 ,y*2 



(3.3) 

T Ict (x) 

Further we differentiate the tail term in the potential (2.10) which arises from source 
contributions interior to the light cone. Apart from the integral 



* i ! ) (?)=9 J drk^["-^)k - = (3.4) 



the gradient $taii,^ = + $m ' contains also local term 



(3.5) 



which is due to time-dependent upper limit of integral in eq.(2.10). Argument of Bessel's 
function W := k oy /-(K ■ K). 

To simplify the tail contribution us much us possible we use the identity 

*> 1 dW (3.6) 



y/-(K ■ K) (K-u) dr 

in the integral (3.4) and perform integration by parts. On rearrangement, we add it to 
the expression (3.5). The term which depends on the end points only annuls fji \ Finally, 
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the gradient of tail term of the potential (2.10) becomes 



T ICt (x) 



= 9 



drk 2 



Ji(W) 
W 



1 + (K ■ a) K _ 



(3.7) 



The particle's position, velocity, and acceleration under the integral sign are evaluated at 
instant r < r rct (x). The invariant quantity 



-(K-u) 



(3.* 



is an affine parameter on the time-like geodesic that links x to z(r); it can be loosely 
interpreted as the time delay between x and z(r) as measured by an observer moving 
with the particle. (For such hypothetic observer the particle is momentarily at rest: its 
four-velocity u a = (1, 0, 0, 0).) 

Because of asymptotic behavior of the first order Bessel's function with argument 
W := k ^-{K -K) the retarded field 



dip rct (x) 
dx>* 

l + (K-a) 
-9 -3 



(3.9) 



T ICt (x) 



I 



9 M 



is finite on the light cone where W = 0. It diverges on the particle's trajectory only. 
Indeed, if the point of emission z(t) approaches the field point x, all the components of 
four-vector K = x — z(t) becomes infinitesimal ones and the distance (3.8) tends to zero. 

The advanced field 



$ adv (x) = 



= -9- 



dip^ix) 
dx>* 

1 + (K ■ a) 



(3.10) 



K,-g 



u. 



(r ad v)' 



+ g 



/ 



drkl 



Ji(W) 
W 



T adv (x) 



1 + (K ■ a) m m 
~2 K ^ ~ ~ 



depends on particle's future history (see Figure 2). Particle's characteristics in the direct 
part are referred to the advanced instant r adv (a;). 



4 Scalar radiation 

The action (2.1) is invariant under infinitesimal transformations (translations and ro- 
tations) which constitute the Poincare group. According to Noether's theorem, these 
symmetry properties yield conservation laws, i.e. those quantities that do not change 
with time. Outgoing radiation removes energy, momentum, and angular momentum from 
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the source. The quantities are defined in standard way [19], as flows of the stress-energy 
tensor density T 

Pldr) = [ da M T s f , (4.1) 



and its torque 

CM = J da a {x»T™ - xTf) , (4.2) 

that flow across a space-like surface £ which intersects a world line at point z(r); d<r M is 
the vectorial surface element on S. While the stress-energy tensor is quite different than 
that in classical electrodynamics: 

dip dp rj^ u i a p dp dp ,2 2 



(see Refs.[18, 20, 21]). To prepare the way for our discussion of a self-force acted on a 
point-like source coupled with massive scalar field, we consider first the relatively simple 
case of a massless scalar field. 



4.1. Massless scalar radiation 

Let parameter ko = 0. The rate of change of the energy- momentum of the retarded field 
is computed by means of the stress-energy tensor (4.3) where field tensor components are 
replaced by the direct field strengths (3.3). We call the part of massless scalar field which 
scales as r _1 the radiative part of field: 

^{x) = -- r {a-k)k,. (4.4) 

(Null vector k^ = K^/r is normalized in such a way that scalar product (k ■ u) = — 1.) 
The part which scales as r -2 does not involve the acceleration: 

= (**-«!.)• ( 4 - 5 ) 

This is the bound part of field. Following Teitelboim [8] , we decompose the stress-energy 
tensor into bound and radiative parts which are separately conserved off the world line. 
To calculate the radiative part it is straightforward to substitute eq.(4.4) into eq.(4.3): 

AnTZ = ^(a-kfk^. (4.6) 



The bound part 



- (-4) t - (-3) 

is the combinations of terms which depend on the retarded distance as r -4 and r~ 3 , 
respectively: 

4irT^ 4) = ^ (k"k v - k»u v - u^k u + u»u v - , (4.8) 

47r7p 3) = ^(a-k) ^k u - W - vTW - tf^ . (4.9) 
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To decompose the angular momentum tensor density 

M^ a = x»T va - x v T m . (4. 10) 

into the bound component and the radiative component, we use the formulae presented 
firstly in pioneer work [22]: 

KZ = ^T^-z'TZ + K^-K'T^, (4.11) 
KZ a = x»TZ-x»TZ + K^ 3) -K»T^ 3) . (4.12) 

We enclose particle's world line by a very thin tube [17] of constant radius r and 
calculate fluxes of the stress-energy tensor (4.1) and its torque (4.2) through this surface. 
An appropriate coordinates are the retarded coordinates [5, Chapt.II] locally given by 

x a = z a (s) +rk a . (4.13) 

Here z(s) G ( is the retarded point associated with a field point x, r is the retarded 
distance between these points, and k a (x) is a null vector field tangent to the congruence 
of null rays that emanate from z(s). 

Globally, the flat space-time M 4 can be thought as a disjoint union of world tubes of 
all possible radii r > 0. A world tube is a disjoint union of spheres of constant radii r 
centered on a world line of the particle. The sphere S(z(s),r) is the intersection of the 
future light cone generated by null rays emanating from point z(s) e ( in all possible 
directions 

C + (z(s)) = {x e M 4 : {x a ~ z a {s)) (x p - z p (s)) = 0, x° - z°{s) > 0} (4.14) 

and the tilted hyperplane 

E(z(s), r) = {x E M 4 : u a (s) (x a - z a (s) - u a (s)r) = 0}. (4.15) 

Points on the sphere S(z(s), r) = C + (z(s)) nH(z(s), r) are distinguished by pair of angles: 
polar angle <d and azimuthal one <p, which specify direction of null vector k (see. [5, 
Chapt.II, Fig.7]). 

A point in Minkowski space can be specified therefore by means of curvilinear coordi- 
nates (s, r, i?, ip). To define a retarded coordinate system, we choose an origin point z(s) 
on particle's world line, future light cone (4.14) with vertex at this point, and tilted hyper- 
plane (4.15) orthogonal to particle's four-velocity u(s). The relation between Minkowski 
coordinates and retarded coordinates is given by eq.(4.13). 

We restrict ourselves to calculation of radiative parts of energy, momentum, and an- 
gular momentum. The technique developed in Ref . [23] can be easily adapted to this task. 
The outward directed surface element d<r M of the cylinder r = const is 

do-^ = r 2 [-Up + (1 + ra^k^] ds A dQ, (4.16) 



11 



where := (a • k) is the component of particle's acceleration in the direction k and 
dQ = sin-^d^Ad^ is the element of a solid angle. The radiative part of energy- momentum 
carried by the massless scalar field is 



T 

J ds J dttr 2 + (1 + rak)^} 9 -{a ■ fc) W 



(4.17) 



dsa 2 (s)u v \s) . 



Integration over angular variables is handled via the relations [23]: 



An 



dttk a = u a , 



(4.18) 



— / dnk a k? = -u a U P + -r] afS , 
An J 3 3 ' ' 

-!- / dnk a kPW = 2u a u^u^ + -u a r]^ + \vPrT + -u 1 ^- 
An J 3 3 3 

The computation of radiated angular momentum is virtually identical to that presented 
above and we do not bother with details. Resulting expression is as follows: 



9_ 
3 



ds { a is] 



z"(s)u u (s) - ^(sK(s) 



(4.19) 



+ u' 1 (s)a v (s)-u v (s)a' 1 (s)}. 



We see that the radiated energy-momentum carried by massless scalar field is equal to one- 
half of the well-known Larmor rate of radiation integrated over the world line. Similarly, 
the radiated angular momentum is equal to the one-half of corresponding quantity in 
classical electrodynamics [22]. This finding is in line with expressions obtained in [13, 15, 
16] for massless scalar self-force. 



4.2. Radiation carried by massive scalar field 

In this Section we decompose the energy-momentum and angular momentum carried by 
massive scalar field into the bound and radiative parts. The bound terms will be absorbed 
by particle's individual characteristics while the radiative terms exert the radiation reac- 
tion. We do not calculate the flows of the massive scalar field across a thin tube around 
a world line of the source. To extract the appropriate finite parts of energy-momentum 
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and angular momentum we apply the scheme developed in Refs.[9, 10]. The scheme sum- 
marizes cumbersome calculations of flows of energy, momentum, and angular momentum 
carried by electromagnetic field of a point-like source arbitrarily moving in flat space-time 
of three dimensions. In 2 + 1 electrodynamics both the electromagnetic potential and the 
electromagnetic field are non-local: they depend not only on the current state of motion 
of the particle, but also on its past (or future) history. The scalar field strengths (3.9) 
and (3.10) behave analogously. 

To find the tail parts of radiated Noether quantities sourced by the interior of light 
cone we deal with the retarded and the advanced fields defined on the world line only: 



^ et = 9 P dr 2 f( Zl ,z 2 ), 

J — oo 

= 9 [ T dr 2 f(z u z 2 ). 



(4.20) 
(4.21) 



The integrand is 



f^{zi,z 2 ) = gk 



Ji[hy/-(Q ■ q)] 

\/-(q ■ q) 



1 + (g • a 2 ) 



U 2„ 



r2 



(4.22) 



where w = k$^J—(q ■ q). Time-like vector g M = z^ — z 2 connects an emission point 
z{t 2 ) G C with a field point z{ji) G (; r 2 = —(q ■ u 2 ). In the forthcoming expressions up 
to the end of the present paper index 1 indicates that the particle's velocity or position is 
referred to the instant r x g] — oo, r] while index 2 says that the particle's characteristics 
are evaluated at instant r 2 before r. 

Let's explain how these fields appeared. In case we do decide to integrate of energy 
and momentum densities over E [9], we should study of interference of outgoing scalar 
waves emitted by different points on particle's world line. This is because massive scalar 
waves propagate not just at a speed of light, but also at all speeds smaller than or equal 
to the speed of light. The world tube is not convenient to study the interference. The 
tilted hyperplane which plays privileged role in the radiation reaction problem in classical 
electrodynamics [8] is not suitable for our purpose too. The reason is that there is no a 
plane which is orthogonal to the particle's 4-velocities at all points on ( before the end 
point z(t) = ( fl S. In Refs.[9, 10] we choose the simplest plane T, t = {x G M 3 : x° = t] 
associated with an unmoving inertial observer. Non-covariant terms arise unavoidable 
due to integration over this surface. To reveal meaningful contribution in radiated energy- 
momentum we apply the criteria which were first formulated in [8, Table 1]: 



• the bound term diverges while the radiative one is finite; 

• the bound component depends on the momentary state of the particle's motion 
while the radiative one is accumulated with time; and 

• the form of the bound terms heavily depends on choosing of an integration surface 
while the radiative terms are invariant. 
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Figure 3: Outgoing scalar waves generated by the portion of the world line that corre- 
sponds to the interval — oo < r 2 < t\ combine within the gray disk in the plane S t = {x G 

M 4 : x° = t}. Corresponding domain of integration is dri dr 2 J^ 1 dR J dQ where 
ki is the radius of smaller disk. After integration over spherical variables the domain 
reduces to double path integral /_ dri f_* dr 2 . It is sketched on the right figure. It is 
obvious that this domain is equivalent to dr 2 dri. 

Since the stress-energy tensor, either electromagnetic or scalar, is quadratic in field 
strengths, we should twice integrate it over ( in order to calculate its flux (4.1) through E t . 
Figures 3 and 4 picture the interference of the disk emanated by fixed point z(ti) G ( with 
radiation generated by portion of the world line from the remote past to the observation 
instant r. Cumbersome calculations performed in [9] can be summarized as a simple 
scheme which can be easily adopted to case of massive scalar field. 

1. Integration of the stress-energy tensor density over E t yields action at a distance 
theory which manipulates with fields evaluated on the world line only. 

2. In case of combination of waves pictured in Fig.3 the surface integration (4.1) of the 
stress-energy tensor contributes in radiated energy-momentum one-half of the work 

P<et = -\£ dnF&fa) (4.23) 

of the retarded tail force (4.20) acting on the charge itself. (The charge "fills" its 
own massive field just as an external one.) 

The surface integration of the angular momentum tensor density (4.2) gives one-half 
of the path integral of the torque of force (4.20): 

MZ = -\ f dn [^fa) - z'FZin)] . (4.24) 
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?2 




Figure 4: Outgoing scalar waves generated by the portion of the world line that corre- 
sponds to the interval T\ < r 2 < r joint together inside the dark disk. The domain of 

integration is dri dr 2 J^ 2 dR J dQ where k\ is the radius of dark disk. After integra- 
tion over spherical coordinates which parameterize E t , the double integral dri dr 2 
survives. It can be replaced by /^ 00 dr 2 J^ 2 dri. 



3. Interference of outgoing scalar waves pictured in Fig.4 takes away from radiated 
energy-momentum one-half of the work 

A** = ~\ f jr.F^r,) (4.25) 

of the advanced tail force (4.21) acting on the charge itself. The surface integration 
of the angular momentum tensor density gives one-half of the path integral of the 
torque of advanced force (4.21): 

= -\ dn WF^in) - z'F^in)} . (4.26) 

4. The the support of double integral dri dr 2 coincides with the support of 
the integral Jl ao dT 2 /I^dTi (see Figs. 3 and 4). Since instants T\ and r 2 label 
different points at the same world line (, one can interchanges the indices "first" 
and "second" in the integrand. Via interchanging of these indices we finally obtain 
J-oo dr i /- 1 o dr 2fl , ( T 2> T i) instead of initial J^dri dr 2 g(T 1 , r 2 ). 

5. The radiative part of energy-momentum carried by massive scalar field is therefore 

PL,» = ~ \ (F^ - F^) (4.27) 

2 rr 



g r r 1 

— dn dr 2 (fnizu z 2 ) - f^(z 2 , Zi)) . 
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Both the radiated angular momentum 



M R 



dri I dr 2 | Zi^/^i, 22) - z^f^Zx, z 2 ) 

I J — oo I. 

[Z2,nfv{z2, Zl) — Z2,uf/j,(z2, Z\) 



(4.28) 



" 2 



and energy- momentum (4.27) exert the radiation reaction. 



Now we evaluate the "upper-limit" behavior of the tail Noether quantities. By this we 
mean the coincidence limit r 2 — > T\ of the of the expressions under the double integrals 
in eqs. (4.27) and (4.28), namely 



fj, 

PtaiLR 



9_ 
' 2 



and 



iW tail,R 



9_ 
2 



f dn /' 

J —oo J —i 



2 JlH 

dr 2 /e 



dr 2 A; - 



1 + (g' fl 2) « _ < , l-(g-Ol) » _ < 



1 + (g • a 2 ) 



i 5 l z i z : 



+ 

Z \Z2 — z l Z2~) + 
2 ~ Z 1 Z 2 ) + 



n 



.29) 



Z^M 2 ^i M 2 



^2 



Z 2 U \ Z 2 U 1 



. (4.30) 



Let T\ be fixed and t\ — r 2 
for our purposes 



A be a small parameter. With a degree of accuracy sufficient 



«2 



A 
A 



(4.31) 



Uj 1 — a^A + 



<- + d 



Substituting these into integrands of the double integrals of eqs. (4.29) and (4.30) and 
passing to the limit A — > yields vanishing expression. Hence the subscript "R" stands 
for "regular" as well as for "radiative" . 

In the specific case of a uniformly moving source g M = u^{t\ — r 2 ) and r\ = r 2 = T\ — r 2 . 
It immediately gives f tl (zi,Z2) = fn(z2,zi) and, therefore, the integrands in eqs. (4.29) 
and (4.30) are identically equal to zero. The local parts of radiation (4.17) and (4.19) 
vanish if = const. As could be expected, nonaccelerating scalar charge does not radiate. 

In the following Section we check the formulae (4.27) and (4.28) via analysis of energy- 
momentum and angular momentum balance equations. Analogous equations yield correct 
equation of motion of radiating charge in conventional electrodynamics [24] as well as in 
flat spacetime of six dimensions [25]. It is reasonable to expect that conservation laws 
result correct equation of motion of point-like source coupled with massive scalar field 
where radiation back reaction is taken into account. 
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5 Balance equations 



The equation of motion of radiating pole of massive scalar field was derived by Harish- 
Chandra [13] in 1946. (An alternative derivation was produced by Havas and Crownfield 
in [26].) Following the method of Dirac [7], Harish- Chandra enclosed the world line of the 
particle by a narrow tube, the radius of which will in the end be made to tend to zero. 
The author calculates the flow of energy and momentum out of the portion of the tube in 
presence of an external field. The condition was imposed that the flow depends only on 
the states at the two ends of the tube (the so-called "inflow theorem", see [27, 28]). After 
integration over the tube along the world line and a limiting procedure, the equation of 
motion was derived. In our notation it looks as follows: 

«*>«? - J « - «?<) - + 9 2 f &< Jj ^<f + 9 2 ^ U f Askl Jl[w) 

= s-r^+^K^,) (5.D 

where mo is an arbitrary constant identified with the mass of the particle and </? ex t is 
the scalar potential of the external field evaluated at the current position of the particle. 
J 2 (u>) is the second order Bessel's function of w — k$^J—(q ■ q). In this Section the 
Harish-Chandra equation will be obtained via analysis of energy-momentum and angular 
momentum balance equations. 

In previous Section we introduce the radiative part = pdir.R + Ptaii.R of energy- 
momentum carried by the field. It defines loss of energy and momentum due to scalar 
radiation. The bound part, p$, is absorbed by particle's 4-momentum so that dressed 
scalar charge would not undergo any additional radiation reaction. Already renormal- 
ized particle's individual four-momentum, say p par t, together with p R constitute the total 
energy-momentum of our composite particle plus field system: P = p pa rt+PR- We suppose 
that the external scalar force g& ex t matches the change of P with time: 

PU+Pk = A- (5-2) 
The overdot indicates differentiation with respect to proper time parameter r. 

Differentiating the tail contribution (4.29) in radiated energy-momentum, we obtain 
the single path integral: 

9 2 [ T r i + (g-a a ) < l-(q-a T ) uf 



w 



(5.3) 



Here index r indicates that particle's position, velocity, or acceleration is referred to the 
observation instant r while index s says that the particle's characteristics are evaluated 
at instant s < r. The expression should be added to the proper time derivative of energy- 
momentum (4.17) carried by a massless scalar field. Substituting the sum p dir +Ptaii f° r 
Pr into eq.(5.2), we obtain the energy-momentum balance equation: 



P^(r) = -^« 2 (r)< + ^ f d^ 
3 2 i 



w 



+ (g • a s) u _ < , l-(g-flr) » 

+ A- (5.4) 
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Our next task is to derive expression which explain how four-momentum of "dressed" 
scalar charge depends on its individual characteristics (velocity, position, mass etc.). 

We do not make any assumptions about the particle structure, its charge distribution 
and its size. We only assume that the particle four-momentum p part is finite. To find out 
the desired expression we analyze conserved quantities corresponding to the invariance 
of the theory under proper homogeneous Lorentz transformations. The total angular 
momentum, say M, consists of particle's angular momentum z A p part and radiative part 
of angular momentum carried by massive scalar field: 

= 4pU(t) - ^(r) + Mf(r). (5.5) 

The radiated angular momentum M£ v = M^ AiT + is determined by eqs. (4.19) 

(direct part) and (4.30) (tail part), respectively. We assume that the torque g{z!^^ xt — 
z^$g Xt ) of the external force matches the change of M with time. Having differentiated 
the angular momentum expression (5.5) and inserting eq.(5.4), we arrive at the equality 

U T A ( Ppalt + 9 4*r+ 9 4 [ T d8k*^±) = 0, (5.6) 

where r T = —(q ■ u(t)). Symbol A denotes the wedge product. 

Apart from usual velocity term, the 4-momentum of dressed scalar charge contains 
also a contribution from field: 

Pnart = mU r ~ TT a r _ "7T / dsfcn . (5.7) 

part 3 2 7.^ w r T v 1 

The local part is the scalar analog of Teitelboim's expression [8] for individual 4-momen- 
tum of a dressed electric charge in classical. The tail term is then nothing but the bound 
part of energy-momentum carried by the massive scalar field: 

KUs = ~\ dri (ifet + FZJ . (5.8) 

The bound part of the field energy-momentum is permanently "attached" to the charge 
and is carried along with it. It is worth noting that the expression under the integral 
sign does not diverge if s — > r. The "local" Coulomb-like infinity is the only divergency 
stemming from the pointness of the source (see Appendix A). 

In the specific case uniform motion = const and argument of Bessel's function 
simplifies: w = k (r — s). Similarly jr T = u^. Since 



r ^Mr-s)} = 

J-oo t-s 



(5.9) 



the field generated by a uniformly moving charge contributes an amount p^ ail s = — 1/ 2g 2 koU tl 
to its energy-momentum. This finding is in line with that of Appendix A where is estab- 
lished that if the particle is permanently at rest, the scalar meson field adds —l/2g 2 k to 
its energy. 
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The expression for the scalar function m(r) is find in Appendix B via analysis of differ- 
ential consequences of conservation laws. We derive that already renormalized dynamical 
mass m depends on particle's evolution before the observation instant r: 



m = m + g 2 I dsfcg ^ ^ - g(p cxt - (5.10) 

J-oc W{T,S) 

The constant mo can be identified with the renormalization constant in action (2.1) which 
governs the dynamics of a point-like charge coupled to massive scalar field. m absorbs 
Coulomb-like divergence stemming from local part of potential (2.10). It is of great 
importance that the dynamical mass, m, will vary with time: the particle will necessarily 
gain or lost its mass as a result of interactions with its own field as well as with the 
external one. The field of a uniformly moving charge contributes an amount g 2 ko to its 
inertial mass. 

To derive the effective equation of motion of radiating charge we replace Pp art in left- 
hand side of eq.(5.4) by differential consequence of eq.(5.7). We apply the formula 

At the end of a straightforward calculations, we obtain 



ma^ + mu^ = (d£ - a*<) + g 2 dsk 2 



3 r r, ■ o u w 



l + (g-O a ) n K 
9 y 



+ g^ xt , (5.12) 



where dynamical mass m(r) is defined by eq.(5.10). The direct part of the self-force 
is one-half of well-known Abraham radiation reaction vector, while the tail one is then 
nothing but the tail part of particle's scalar field strengths (3.9) acting upon itself. 

Now we compare this effective equation of motion with the Harish- Chandra equation 
(5.1). The latter can be simplified substantially. Having used the recurrent relation 

J 2 H = J J™1- (5.13) 
w aw 

between Bessel functions of order two and of order one, after integration by parts we 
obtain 



9 1 f dsk^ q » = ^ 2 < - 9 2 f dskr-^- i^£V - ^ • (5.14) 



W 



l + (g-a«) „ _ < 

.2 " 



rt r 



We also collect all the total time derivatives involved in Harish- Chandra equation (5.1). 
The term m(r)u^ arises under the time derivative operator, where time-dependent func- 
tion m(r) is then nothing but the dynamical mass (5.10) of the particle. On rearrange- 
ment, the Harish- Chandra equation of motion (5.1) coincides with the equation (5.12) 
which is obtained via analysis of balance equations. It is in favor of the renormalization 
scheme for tail theories developed in the present paper. 
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To clear physical sense of the effective equation of motion (5.12) we move the velocity 
term mu^ to the right-hand side of this equation: 



(5.15) 



According to [14], the scalar potential produces the Minkowski force 

= g (»r + «) ^ 

which is orthogonal to the particle's 4-velocity. The self-force 



J —c 



dsk^ 
w 



1 + (g • a «) _ ^^M) _ < + K ■ 



(5.16) 



(5.17) 



is constructed analogously from the tail part of gradient (3.9) of particle's potential (2.10) 
supported on the world line (. Indeed, since the massive scalar field propagates at all 
speeds smaller than the speed of light, the charge may "fill" its own field, which will act 
on it just like an external field. The own field contributes also to particle's inertial mass 
m(r) defined by eq.(5.10). 



6 Conclusions 

In the present paper, we adopt the Dirac scheme of decomposition of the retarded Green's 
function into symmetric (singular) and radiative (regular) parts to functions supported 
within light cones. The regularization scheme summarizes a scrupulous analysis of energy- 
momentum and angular momentum balance equations in 2 + 1 electrodynamics [9, 10]. It 
differs from the approach developed by Detweiler and Whiting [6] on two "extra" entities: 
additional instant T\ before the instant of observation and extra integration of the retarded 
and the advanced tail forces over particle's path. So, the retarded tail force depends on 
the particle's past history before T\ < t (see Figure 5). Its advanced counterpart is 
generated by portion of the world line that corresponds to the interval [ti,t]. The tail 
part of radiated energy-momentum is one-half of the work done by the retarded force 
minus one-half of the work done by the advanced force, taken with opposite sign. This 
part of radiation detaches the point source and leads an independent existence. 

The support of the advanced force is a portion of the world line that corresponds to 
finite time interval. It is presented implicitly in the Detweiler and Whiting construction 
where the tail term (1.6) is introduced to achieve appropriate retarded causality. It 
"evaporates" in the coincidence limit x — > z(t) at which its support [r rct (x) , r adv (x)} 
shrinks to zero. 

The one-half sum of the retarded and the advanced works is the bound part of tail 
energy-momentum which is permanently attached to the charge and is carried along with 
it. It modifies particle's individual characteristics (its momentum and its inertial mass). 
A point source together with surrounded "cloud" constitute dressed charged particle. 
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T 2 



Figure 5: We call "retarded" the force (4.20) with integration over the portion of the 
world line before T\. We call "advanced" the force (4.21) with integration over the portion 
of the world line after T\. For an observer placed at point z{r\) G ( the radiative (4.27) 
and the bound (5.8) parts of massive scalar field momentum looks as the combination of 
incoming and outgoing radiation, and yet the retarded causality is not violated. We still 
consider the interference of outgoing waves presented at the observation instant r. The 
scalar field carries information about the charge's past. 

Since the properties of the retarded and the advanced solutions of wave equation, the 
one-half sum is singular while one-half difference is regular at the location of the particle. 

The bound and the radiative angular momentum carried by charge's field are simply 
torques of the above combinations of the retarded and the advanced tail forces. 

Together with contributions from the direct part of Green's function, the tail terms 
constitute Noether quantities removing by outgoing waves from the dressed source. Chan- 
ges in individual momentum and angular momentum of dressed charged particle compen- 
sate losses of energy, momentum, and angular momentum due to radiation. (Influence 
of an external device can be modelled easily.) Analysis of balance equations yields the 
Harish-Chandra equation of motion of radiating scalar pole [13]. This equation includes 
the effect of particle's own field as well as the influence of an external force. 

Energy-momentum and angular momentum balance equations for radiating scalar pole 
constitute system of ten linear algebraic equations in variables Pp art (T) and their first time 
derivatives Pp ar t( r ) as functions of particle's individual characteristics (velocity, accelera- 
tion, charge etc.). The system is degenerate, so that solution for particle's 4-momentum 
includes arbitrary scalar function, m(r), which can be identified with the dynamical mass 
of the particle. Besides renormalization constant, the mass includes contributions from 
particle's own field as well as from an external field. 

This is a special feature of the self force problem for a scalar charge. Indeed, the 
time- varying mass arises also in the radiation reaction for a pointlike particle coupled to a 
massless scalar field on a curved background [4]. The phenomenon of mass loss by scalar 
charge is studied in [29, 30]. Similar phenomenon occurs in the theory which describe a 
point-like charge coupled with massless scalar field in flat spacetime of three dimensions 
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[31]. The charge loses its mass through the emission of monopole radiation. 
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Appendix A Energy-momentum of the scalar massive field of uniformly mov- 
ing source 

In this Appendix we calculate the energy-momentum (4.1) carried by the scalar massive 
field due to static charge g. The stress-energy tensor T is given by eq.(4.3) where tp is 
Yukawa potential (2.13). 

It is convenient to choose the simplest plane E t = {x G M 4 : x° = t} associated with 
unmoving observer. We start with the spherical coordinates 

x° = s + r, x % = rn l (6.1) 

where ri 1 = (cos (j> sin 6, sin sin 9, cos 9) and s is the parameter of evolution. To adopt 
them to the integration surface S t we replace the radius r by the expression t — s. On 
rearrangement, the final coordinate transformation (x° . Ob . Ob . Ob 3 ) !->■ (t,s,(f>,9) looks as 
follows: 



x 



t, x l = (t - s)n l 



The surface element is given by 



(6.2) 
(6.3) 

where dfl = sin #d#d0 is an element of solid angle. 

After trivial calculation one can derive the only non-trivial component of energy- 
momentum (4.1) is 



da = (t - sfdsdtt 



Pi 



dn- 

2 



V ( — 



1 2 2 



(6.4) 



2 



k exp[-2k (t - s)] + 



exp[— 2k (t — s)] 



t-s 



s-tt 



9" 9' 



= lim 

e^o2e 2 

where e is positively valued small parameter. 

Having performed Poincare transformation, the combination of translation and Lorentz 
transformation, we find the energy-momentum carried by massive scalar field of uniformly 
moving charge: 

p£ = hmJ-^-^U ^- (6.5) 

The divergent Coulomb-like term is absorbed by the "bare" mass m involved in action 
integral (2.1) while the finite term contributes to the particle's individual 4-momentum 
(5.7). 
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Appendix B Derivation of the renormalized mass of scalar charge 



We follow the scheme elaborated within analysis of 2 + 1 electrodynamics where the ex- 
pression for the scalar function m(r) is derived via analysis of differential consequences 
of conservation laws. In hypothetical Minkowski space of three dimensions already renor- 
malized mass of charged particle depends on particle's evolution before the observation 
instant r. 

The scalar product of particle 4- velocity on the time derivative of particle 4-momentum 
(5.4) is as follows: 



/ . \ g 2 2 g 2 
(ppart-ur) = y^ + y 

+ #(w$ cxt ). 



dsk; 



JAw) 



w 



1 + (g • as) ( \ K • U r) , (g • gr) 



(6.1) 



Since (u • a) — 0, the scalar product of particle acceleration on the particle 4-momentum 
(5.7) does not contain the scalar function m: 



, \ g 2 2 g 2 

(Ppart • Or) = - — 



f 



dsk 



2 J 1 (w) (q ■ Or) 



w 



Summing up (6.1) and (6.2) we obtain 



d 

d7 



(Ppart ' U T ) 



g_ 

2 



dsk? 



JiH d 
w ds 



(q ■ u 7 



(6.2) 



(6.3) 



We rewrite the expression under the integral sign as the following combination of partial 
derivatives in time variables: 



k? 



J\{w) g f (q ■ u T ) 

w ds 



d_ 

ds 



2 Ji(w) (q ■ U T ) 



w 



d_ 



2 Ji(w) 



w 



(6.4) 



At this point we suppose that the external field is the gradient of an external scalar 
potential, say ip cxt (x). If the field is referred to the point z(t) where the charge is located, 
the scalar product (u ■ $ cx t) becomes the total time derivative: 



(u ■ $ 



dz a dip 



cxt V 



cxt j 



dr dz a 
dv?cxt 



(6.5) 



dr ' 

These circumstances allow us to integrate the expression (6.3) over r: 



(Ppart -u T ) = -mo + 



dri 
dr 2 



dr 2 



dn 



_d_ 

dr 2 



Ji(w) (q ■ Mi) 



w 



r 2 



T2 



_d_ 

= -m 



k 



2 Ji{w) 



w 



+ g^ext 



I dsk 2 



Ji(w(t,s)) 
w(t, s) 



+ gv 



cxt- 



(6.6) 
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The external potential is referred to the particle's position z(t). 

Alternatively, the scalar product of 4-momentum (5.7) and 4-velocity is as follows: 

GW ■u T ) = -m+ 9 ^ f dse Jl[w . {T '° )] . (6.7) 
2 w(t,s) 

Having compared these expressions we obtain: 

m = mo + g 2 f T dskt Jl[w{T > S)] - g^ x Mr)]- (6.8) 

J-oo W(T,S) 

The constant to can be identified with the renormalization constant in action (2.1) which 
governs the dynamics of a point-like charge coupled to massive scalar field, too absorbs 
Coulomb-like divergence (6.5) stemming from direct part of potential (2.10). It is of great 
importance that the dynamical mass, to, will vary with time: the particle will necessarily 
gain or lost its mass as a result of interactions with its own field as well as with the external 
one. Since eq.(5.9), the field of a uniformly moving charge contributes an amount g 2 k to 
its inertial mass. 
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